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Recent years have seen increased interest in complexified Bohmian mechanical trajectory calcu- 
lations for quantum systems, both as a pedagogical and computational tool. In the latter context, 
it is essential that trajectories satisfy probability conservation, to ensure they are always guided 
to where they are most needed. In this paper, probability conservation for complexified Bohmian 
trajectories is considered. The analysis relies on time-reversal symmetry considerations, leading to 
a generalized expression for the conjugation of wavefunctions of complexified variables. This in turn 
enables meaningful discussion of complexified flux continuity, which turns out not to be satisfied in 
general, though a related property is found to be true. The main conclusion, though, is that even 
under a weak interpretation, probability is not conserved along complex Bohmian trajectories. 



I. INTRODUCTION 



Trajectory interpretations of quantum mechanics have 
been of interest since the earliest days of the quantum 
theory. Indeed, they even predate the Schrodinger equa- 
tion itself — as one finds, e.g., by considering the Bohr- 
Sommcrfeld quantization rule. Although the latter was 
discovered to be an incorrect description of quantum the- 
ory, it survives today in the form of the Jeffrey- Wentzel- 
Kramers-Brillouin (JWKB) approximation, or more gen- 
erally, semiclassical mechanics In this approach, a 
time-evolving quantum pure state is treated as a statisti- 
cal ensemble of classical trajectories that "carry" approx- 
imate quantum information, i.e. complex amplitudes. 
There are many reasons, both pedagogical and practical, 
why semiclassical and even classical trajectory methods 
may be regarded as beneficial. Surely though, one of 
these must be the fact that the differential probability, 
p(x)dx, is conserved along any given trajectory — a well- 
known property of classical statistical ensembles. This 
alone ensures that the ensemble trajectories travel to 
where they are "most needed," i.e. to where the probabil- 
ity density is largest — a consideration that is especially 
important for localized wavepacket propagation in the 
limit of large system dimensionality. 

Starting with Madelung in the same year as the 
Schrodinger equation itself Q (based on matter wave 
ideas of deBroglie) and evolving into a full-fledged in- 
terpretation of quantum theory with D. Bohm and 
coworkers in the early 1950's [H, H|, an exact trajec- 
tory formulation of quantum mechanics has also been 
developed. Over the ensuing decades, the resultant 
"Bohmian mechanical" trajectory ensembles have been 
relied upon to provide interpretational or "analytical" 
insight into solved time-dependent quantum wavepacket 
propagation problems, such as the fundamental double- 
slit experiment [1, HI 0, fi Q- More recently, innova- 
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tions spearheaded by members of the chemical physics 
community have led to the use of quantum trajectory 
methods (QTMs) as a "synthetic" tool, i.e. to solve 
the time-dependent Schrodinger equation (TDSE) itself 

Hamuli. 

Though it is of great interest to compare and contrast 
the behavior of quantum trajectories with their classical 
counterparts, we shall do so here only as it relates to 
the present goals, as a detailed discussion would take us 
too far afield. Most important in the present context is 
the fact that Bohmian quantum trajectories also satisfy 
probability conservation [Eq. (|9|)] — an extremely benefi- 
cial property for the "synthetic" application of QTMs, 
and again, a chief reason for their utility. On the other 
hand, standard Bohmian mechanical trajectories — which 
we shall henceforth refer to as "real-valued Bohmian tra- 
jectories" for reasons that will become clear — suffer from 
certain apparent drawbacks as well, some of which can 
cause severe numerical difficulties for the synthetic ap- 
proach [Toj |. In particular, for non-degenerate station- 
ary states, all quantum trajectories are stationary fixed 
points — in stark contrast to the corresponding classical 
trajectory orbits. 

To circumvent the above problems, one approach is 
to follow the semiclassical prescription of adopting a 
multipolar expansion of the wavefunction, [131 Il4j . 
This leads to correspondence between individual (real- 
valued) quantum and classical trajectories in the classical 
limit. A second approach, the focus of the present pa- 
per, involves a different generalization of Bohm's original 
prescription — i.e., allowing the coordinates, x, and tra- 
jectory velocities, v, to take on complex values. The re- 
sultant "complex-valued Bohmian trajectories" offer cer- 
tain advantages; for instance, they are known not to 
be fixed-points, in general, for nondegenerate station- 
ary states, so that it is possible to achieve nontrivial 
trajectory dynamics in this context. Although complex- 
valued Bohmian mechanics may still be in its infancy, 
interest has grow n tremendously in the last few years 
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the 1980's with a paper by Leacock and Padgett [l5| and 
a less well known (and very brief) article by Tourcnnc 
[T3 |. More recent authors have explored the complex 
Bohmian approach both for time-independent (station- 
ary) and time-dependent (wavepacket propagation) prob- 
lems, in both the analytical and synthetic contexts (l7L 
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Remarkably, the all-important issue of trajectory prob- 
ability conservation does not yet appear to have been ad- 
dressed in the literature. Yet as stated previously, this 
is an essential requirement if the complex-valued syn- 
thetic TDSE methods currently under development are 
to have general utility for systems larger than two or 
three dimensions, d. These methods already require spe- 
cial root-finding procedures to single out the subset of 
complex trajectories that happen to arrive on the d- 
dimensional real "axis" a; at a desired final time t / — since 
at other times, this trajectory ensemble is described by 
some non-trivial <i-dimensional manifold embedded in the 
2e?-dimcnsional complex coordinate space. Alternatively, 
one could in principle avoid the root search by simulta- 
neously propagating trajectories over the entire complex 
coordinate space, but only at the expense of doubling 
the dimensionality of the trajectory ensemble space, and 
thus enormously increasing the required number of tra- 
jectories to compute. Such an approach would be highly 
redundant from an information theory perspective, but 
perhaps feasible much like coherent state representa- 
tions [ij. 

Though daunting, it may bepossible to overcome the 
above "d-doubling problem" [3I| for large systems, but 
only if complex quantum trajectories turn out to satisfy 
probability conservation, at least approximately. This 
paper thus considers the issue of probability conserva- 
tion for complexified Bohmian trajectories. The analysis 
ultimately relies on time-reversal symmetry considera- 
tions, leading to a generalized definition of the conju- 
gation operation for functions of complexified variables, 
which in turn leads naturally to a discussion of complex- 
ified flux continuity and probability conservation. An- 
alytical properties of complex trajectory dynamics are 
then considered, and finally, several d = 1 examples are 
discussed. 



II. BACKGROUND 

Real-valued Bohmian mechanics begins with the 
Madelung-Bohm ansatz for the wavefunction 0, d, 0| > 



V(x,t) =R(x,t)e lS{x < t)/n . 



(1) 



The decomposition above is essentially unique, by virtue 
of the fact that R and S are taken to be real-valued, 
and R > 0. By substituting Eq. Q into the TDSE, 
and gathering real and imaginary terms separately, one 
obtains the following two real-valued equations: 

dR 



dt 



-1 
2m 



(2R'S' + RS") 



(2) 



dS_ 

dt 



2m ' 2m R 



(3) 



where the primes denote spatial differentiation. Equa- 
tion ([3]) is the "quantum Hamilton- Jacobi equation" [3[, 
which imparts a classical-field-theory-like [33j interpre- 
tation to quantum wavepacket propagation, provided 
that: (a) S'(x,t) is interpreted as trajectory momen- 
tum; (b) an additional "quantum potential," Q(x,t) = 
— (H 2 /2m) (R"/R), is added to the true potential, V(x), 
to determine the trajectory dynamics. 

Equation ([2]) is the flux continuity equation, which un- 
der the above interpretations, is identical to trajectory 
probability conservation. In particular, Eq. can be 
rewritten as 



dp(x,t) 
dt 
p(x,t) 

j(x,t) 
v{x, t) 



-j'(x,t), where (4) 

R 2 (x,t) is probability density, (5) 
p(x,t)v(x,t) is probability flux, and (6) 
S'(x,t)/m is trajectory velocity. (7) 



Note that the flux j(x,t) above accords with the usual 
quantum definition, i.e., 
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2im 



(**4f'-*'*^) 



(8) 



Note also that the flux is independent of potential energy 
as is reasonable, i.e. V(x) enters only into the dynam- 
ical equation for trajectory evolution, Eq. ([3]). In any 
case, Eq. (01 implies conservation of differential proba- 
bility along a trajectory, i.e. 



d[p{x, t)dx] 
Jt 



= 0, 



(9) 



where d/dt refers to the total (hydrodynamic) time 
derivative. For simplicity, the above equations have been 
written as if the dimensionality is d = I, but they are 
meant to refer also to the multidimensional case, as the 
d > 1 generalizations are straightforward. 

In complex Bohmian mechanics, an altogether different 
ansatz is employed: 



$(x,t) = e 



_ iS(x,t)/H 



(10) 



In effect, the real-valued amplitude R of the Eq. (Q]) de- 
composition is "subsumed" into the exponent to form 
the imaginary part of the now- complex- valued action, 
S. Substitution into the TDSE now yields the single, 
complex-valued equation, 



dS_ 
~dt 



2m 



V 



2m 



-S" 



(11) 



which can be interpreted as a complex quantum Hamil- 
ton Jacobi equation. Equation (jf I p is the starting point 
for the quantum Hamilton- Jacobi formalism, regarded 
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as one of the nine fundamental formulations of quantum 
mechanics, pi ITl l34j 

Unlike the real- valued Eq. |3]), Eq. (fTTj) contains all 
of the information present in 'J, leading some authors 
to conclude that the complex version is the more funda- 
mental [HI]. On the other hand, in going from real- to 
complex- valued formulations we have apparently lost the 
flux continuity relation altogether! Thus, if there is in- 
deed a probability conservation property for complex tra- 
jectories, it is neither manifest in, nor independent from, 
Eq. (fTTj) . Note that in certain contexts, it is possible to 
extract a complex energy-momentum conservation rela- 
tion from Eq. (|TT|) pH35|. However, this is not directly 
useful for synthetic applications, vis-a-vis the guidance 
of trajectories to high-probability regions where they are 
needed most. 



III. COMPLEXIFICATION: TIME REVERSAL 
SYMMETRY AND GENERALIZED COMPLEX 
CONJUGATION 

To make progress with regard to probability conserva- 
tion, it seems clear that Eq. (flT|) must be written in a 
different, but mathematically equivalent, form. In com- 
paring Eq. (TJJ) to Eq. (fTO)) . which emphasizes the action S 
over amplitude R, it seems clear that the opposite proce- 
dure should be applied to emphasize the latter — i.e. the 
real-valued S of Eq. (TJJ should be "pulled down" from 
the exponent to form a complex- valued amplitude, R. 
But this would simply yield W itself as the appropriate 
quantity to work with. Thus, the TDSE in its standard 
form (i.e. in terms of should be a reasonable starting 
point for analyzing flux continuity — albeit a complexified 
version, with x S 1Z replaced with zeC. Complexifica- 
tion of the TDSE per se offers no inherent difficulty pro- 
vided that V(x) and the initial wavepacket, ^(a;,i = 0), 
are both analytic functions, as everything can then be 
uniquely lifted from the real x axis to the complex z 
plane in accord with the usual rules of analytic continua- 
tion [36|. On the other hand, establishing a complexified 
probability density, does pose a bit of a problem, 

as the traditional definition p(x) = ^{x)^{x) in terms 
of the conjugate function <5>(x) — ^*(x) — [i'(x)]* is in 
general not analytic when x is replaced with z. 

Two straightforward candidates for a complexified p[z) 
quantity are presented below: 

1. Define p(z) = **(z)*(z). 

2. Define p(z) as the analytic continuation of p{x). 

Option 1. has the presumed advantage that p(z) is pos- 
itive and real- valued everywhere in the complex plane, 
but suffers from the severe drawback that p(z) is not an 
analytic function. This approach is being considered by 
other authors (32j. Option 2. offers analyticity, but yields 
complex-valued probability densities off of the real axis, 
and bypasses the more fundamental issue of non-analytic 



conjugate wavefunctions, which must still be resolved. 
Analyticity is an enormous advantage, for it means that 
familiar expressions such as f_°° p(z)dz — 1 represent 
true contour integrations with path-independent mean- 
ing. Thus, integration contours may be deformed away 
from the real axis — which is especially important for syn- 
thetic TDSE applications, given that the contour is es- 
sentially the (time-evolving) trajectory ensemble mani- 
fold. In any case, complex probability values are likely 
unavoidable in complexified space, for even if the den- 
sity p(z) is real-valued, the differential probability itself, 
p(z)dz, need not be. 

We will resolve the matter by directly addressing the 
more general and fundamental issue of wavefunction con- 
jugation on complexified spaces. One approach to this 
problem is to invoke charge/parity/time-reversal (CPT) 
symmetry [37j — an idea that was introduced previously 
in the specific context of non-Hermitian Hamiltonian op- 
erators [H, [H, H^] . Similar ideas can be applied in the 
present case of complexified Hermitian Hamiltonians, al- 
though for general potentials, only time-reversal symme- 
try is relevant. Let ^(z, t) be a solution of the complex- 
ified TDSE (z e C, but t € K): 



ih 



dt 



h z <9 2 4<(z,i) 
2m dz 2 



(12) 



Complex conjugating both sides and applying explicit 
time-reversal, i.e. t — > —t, yields 



. ft a [*(*,-*)]' 
dt 



2m 



d 2 V(z,-t) 



dz 2 

[V(z)]*[*(z,-t)} 



(13) 



The above is the usual means of deriving the effect of 
the antiunitary time-reversal operator on a wavefunction, 
for real- valued coordinates [37|. For the complexified 
case, however, it is clear that Eq. (|12p is not equivalent 
to Eq. [fig]), because [V (.?)]* = V*(z) is not equivalent 
to V(z) off of the real axis. Instead, we must introduce 
the additional and final step of replacing the coordinate 
z with its complex conjugate, z*. From the Schwartz re- 
flection principle (3q |. V*(z*) = V(z), because V{x) is 
real- valued. Thus, 



dV*(z*,-t) _ h 2 dH*(z*,-t) 



dt 



2m 



dz 2 



+ V(z)**(z*,-t), 
(14) 

so &*(z*, —t) is also a solution of the TDSE, the proper 
manifestation of time-reversal symmetry on the complex- 
ified space. Note that the final step above is also consis- 
tent with the approach of Hiiber, Heller, and Littlejohn 
pfj| |. which may be regarded as a semiclassical approx- 
imation to complexified Bohmian Mechanics. In fact, 
there is a very close connection between the present ap- 
proach and complex semiclassical mechanics, with the 
latter emerging as the first-order truncation of an infinite 
series expansion of the former, as elucidated by Tannor 
and coworkers [H HI Hi, IH . 
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Following Bender's approach [38l 1391] . the correspond- 
ing conjugation operation is therefore 

(15) 



f( Z ) = r(z*)=g( Z ). 

Note that along the real axis, i.e. z = x, Eq. (fl5|) is 
equivalent to the standard conjugation operation. Most 
importantly however, we find that /*(z*) is an analytic 
function in the original variable z, in that it satisfies the 
Cauchy-Reimann conditions [36| — unlike, say, f(z*) or 
/*(z). Equation (| 1 5[) also satisfies the mathematical def- 
inition of a conjugate linear map [36j], unlike f(z*) or 
f*(z). Note that /*(z*), is not equivalent to f(z) in 
general. We will therefore sometimes refer to f*(z*) as 
"3(2)," to emphasize both its distinctness from f(z) and 
also its analyticity. The latter property can be demon- 
strated explicitly from a Taylor expansion: 



E 

k=Q 



C k z k 



(16) 



fe=0 



Having shown that f(z) is analytic, it follows trivially 
that the complexified inner product integration, 

/+00 
g(z)h(z)dz, (17) 
-00 

is contour path-independent, as the integrand itself is an- 
alytic [assuming analytic h(z)]. In particular, this implies 
that 



p(z) = ^(z)^{z) 



(18) 



is analytic, and moreover, is equivalent to option 2. de- 
scribed above. 

Note that both /(z) = V(z) and f(z) = p(z) share 
the special property that f(z) — g(z) = f(z). This is 
nothing but the Schwartz reflection principle, and is true 
whenever f(z) is analytic and f(x) is real. We find it con- 
venient to refer to such a function as a REAL function — 
even though it is understood that f{z) is not real-valued 
off of the real axis. Similar ly, an IMAG analytic function 
f(z) is defined such that f(z) — g{z) — —f(z), and is 
pure imaginary-valued along the real axis, but not neces- 
sarily elsewhere. Any analytic function /(z) can be de- 
composed into a sum of REAL and IMAG parts, which 
we demonstrate via explicit construction: 

f(z) = /re(z) + /im(z), where (19) 

/(*) + /*(*•)" 



/re (2) = 
fm(z) = 



/(*)-/*(**) 



Note that /re(z) and fwi(z) are themselves analytic — 
unlike say, Re[/(z)] and Im[/(z)] (although /re(^) = 
Re[/(ir)], etc.) Note also that the analytic derivative 
f(z) = df{z)/dz of a REAL(IMAG) function /(z) is 
also REAL(IMAG). Finally, the product of a REAL and 
REAL (or IMAG and IMAG) pair of functions is REAL, 
whereas the product of REAL and IMAG functions is 
IMAG. 



IV. FLUX CONTINUITY AND PROBABILITY 
CONSERVATION 

Our next goal is to define a complexified flux, from 
which to derive a corresponding flux continuity rela- 
tion. We start with the velocity field y(z, t), which from 
Eqs. © and (JTUJ) is given by [H [H El, ill 



v(z,t) 



ih$'(z,t) 
m *(z,t) ' 



(20) 



In general, v(z, t) is neither REAL nor IMAG. Along the 
real axis, vre(x, t) is the "flow velocity" (i.e. the stan- 
dard velocity field of real- valued Bohmian mechanics, 
closely related to hydrodynamics), whereas vim.(x, i) is 
known as the "Einstein osmotic velocity," associated with 
wavepacket spreading, or diffusion in stochastic quan- 
tum mechanics [l(J EE EH EH- However, Viu(x,t) is 
not used in the stochastic context to generate trajecto- 
ries per se, or otherwise venture off of the real coordi- 
nate axis. Hirschfelder did use the imaginary velocity 
to create real-valued trajectories or "streamlines," but 
only along the real axis, by replacing t with it [44l| . In 
complexified Bohmian Mechanics, the imaginary velocity 
is directly responsible for transporting the complex tra- 
jectories off of the real axis. Of key significance for the 
present approach is that both velocity components can be 
given significance off of the real axis, due to the Eq. (I19p 
decomposition, and the fact that i>re(z,£) and viM(z,t) 
are analytic functions, provided v{z, t) is (at least locally) 
analytic. Note that global analyticity of ^(z, t) does not 
necessarily imply the same property for v(z, t) (Sec. FV)) . 

The complexified flux is naturally defined from 
Eqs. (fT8")l and (|20p. and the analytic continuation of 
Eq. ©, as 



j(z,t) = v(z,t)p(z,t) 



if) 



y*(z*)V'(z), (21) 



which based on the multiplication rules given in Sec. IIII[ 
is again neither REAL nor IMAG. Note that the conjuga- 
tion operation of Eq. (fT5j) commutes with spatial differen- 
tiation, i.e. [df(z*)/dz*]* = g'(z), where g(z) = f*(z*). 
Thus, the meaning of expressions such as the following 
are unambiguous: 

JreOM) = v KE (z,t)p(z,t) 

= -— hr (z*)*'(z) - <H(zW*(z*)] (22) 
2m 

ji M (z,i) = vi M (z,t)p(z,t) 

= -— hr (z*)*'(z) + <H(zW*(z*)] (23) 
2m 

Note that along the real axis, Jre(^,^) is equivalent to 
the usual real- valued quantum flux of Eq. ([H]). The imag- 
inary flux, j'im(z, t), in contrast, does not appear to have 
been considered previously in the literature, even when 
restricted to the real coordinate axis. Off of the real axis, 
j and v point in different directions. 
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A flux continuity relation should presumably involve 
the divergence of the flux. From Eqs. (|22[) and (|23[) the 
REAL and IMAG components are found to be: 

Jre(M) = -*(«)*"*(z*)](24) 
2m 

+2tf'*(z*)*'(z)] (25) 

To relate the above to the time derivative of p(z,t) as 
defined in Eq. ([TBI) requires the time-derivative of ^(z, t). 
By replacing t — > — < in Eq. HU), this is found to be 



/7i 



0*(z,t) 



0* 



+— *(*,*) -V(*)*(*,f) (26) 

zm 



Multiplying Eq. ([12]) by -(i/h) i f>(z,t), and adding to 
— (i/?i)\E , (z, t) times Eq. (l2l?l) then yields the complexi- 
fied flux continuity relation 



dp(z,t) 
dt 



(27) 



The above procedure is similar to that used to derive 
the real-valued flux continuity equation; however, it only 
works here by virtue of the Schwartz reflection principle. 
Indeed, a continuity relation based on any definition of 
^(z, t) other than Eq. (|15[) would result in a flux quantity 
that depends on the potential energy — a highly unphysical 
scenario that we reject out of hand. 

On the other hand, the continuity relation of Eq. (12"T1) 
involves only the REAL component of flux, rather than 
j(z, t) itself. This makes perfect sense when one considers 
that p itself is REAL, and therefore its time derivative 
must also be REAL. However, this has nontrivial ram- 
ifications vis-a-vis probability conservation for complex 
trajectories generated from v{z, t) rather than from say, 
UrbC 2 ) t)- We will have more to say on this topic in a mo- 
ment, but first we point out some additional noteworthy 
aspects of Eq. (|27|) . In particular, 



dp(z,t) 
dz 



dz = 0, 



(28) 



i.e. the total probability conservation property, which 
has a contour path-independent meaning; it is in any 
event clear that integration of Eq. (|24|) along the real 
axis is zero. In the specific case of nondegenerate sta- 
tionary states, one can further show that dp(z, t)/dt = 
everywhere. Thus, even though p has a nontrivial com- 
plex phase off of the real axis, this phase does not evolve 
over time, so that p{z, t) is truly stationary everywhere. 
Note that without loss of generality, ^(z,t) = ip(z)e~ tuJt 
with ip( z ) REAL in this case, implying that 



p(z,t) = [</>(z)] 2 



(29) 



Since the continuity relation in Eq. Q27[) is given in 
terms of j're(z, t) = vre( z , t)p(z, t), it immediately fol- 
lows via analytic continuation that d[p(z,t)dz]/dt = 0, 



provided the complex trajectories are obtained from 
wre(z, t) rather than v(z, t). Along the real axis, the v-re 
trajectories are just the standard real-valued Bohmian 
trajectories, so the orbit is the real axis itself. Off of the 
real axis, the ure(<£, i) are quite nontrivial, and might 
in principle be considered for dynamical purposes, with 
a ready-built trajectory probability conservation prop- 
erty. There are at least two drawbacks to this arrange- 
ment however: (1) vre trajectories that start off of the 
real axis never intersect the real axis; (2) nondegenerate 
stationary states still have t>RE(z>i) = everywhere in 
the complex plane. In any event, all of the previous lit- 
erature on complexified Bohmian mechanics uses v(z,t) 
rather than Ure(z> t) trajectories. Note that for nonde- 
generate stationary states, v(z,t) = viu(z,t), so that all 
of the trajectory dynamics is due to the imaginary ve- 
locity in this case. Again, because viM{z,t) is not pure 
imaginary-valued, this does not imply trivial "vertical" 
(parallel to imaginary axis) trajectory orbits in the com- 
plex plane. In general, stationary state orbits are vertical 
only where they intersect the real axis, about which the 
orbits display reflection symmetry; otherwise, they are 
quite arbitrary, often recrossing the real axis again at a 
different point (Sec. [VTJ) . The ground state of the har- 
monic oscillator system, for instance, is characterized by 
concentric circular orbits [n| EE 0, H(| • 

For the more general case of nonstationary wavepacket 
propagation, it is easily shown that d[p(z, t)dz]/dt ^ in 
general, for trajectories obtained from v(z,t). Thus, in a 
literal sense, probability conservation along trajectories is 
not satisfied. On the other hand, it might be argued that 
such a strong form of probability conservation should not 
be required. In particular, this relates to the fundamental 
question of how physical observables are to be interpreted 
in complexified space, and how these relate to real mea- 
surements. The impression one derives from the litera- 
ture is that observables are to have physical meaning only 
when evaluated along the real axis. On the other hand, a 
theory such as that presented here, which allows all quan- 
tities to be analytically continued, may shed additional 
light into this question, in that knowledge of said quanti- 
ties along an essentially arbitrary d-dimensional manifold 
implies knowledge along the real axis. 

Returning to the issue of probability conservation, if we 
presume that p(z,t)dz has no direct physical meaning off 
of the real axis, then it may not be appropriate to impose 
Eq. ([9]) throughout the entire complex plane. Instead, 
let us imagine a particular v(z,t) trajectory, z(t), which 
crosses the real axis at time ti , and again at a later time 
tf, i.e. z(ti) = Xi and z(tf) — Xf (Fig. [T]). We stipulate 
on physical grounds that it should be sufficient to satisfy 
the following weak probability conservation condition: 



p(xi,ti)dxi = p(x f ,tf)dxf 



(30) 



However, even Eq. (|30|) turns out not to be satisfied in 
general — even under the most favorable situation where 
is a non-degenerate stationary state and v(z) is ana- 
lytic. 
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V. ANALYTIC PROPERTIES AND 
TRAJECTORY DYNAMICS 

To demonstrate that Eq. (f3T)|) is not satisfied, it is help- 
ful to introduce the map f lyti ' tf \z i ) — Zf, corresponding 
to the trajectory that connects the initial point [z{, U) to 
the final point [zf, tf) (Fig. [l}. This in turn requires the 
trajectory guidance equation, 

dx ih *'(z,t) 

— =v(x,t) = —, r- (31) 

dt y ' m $(z,t) v ; 

Regarded as a function of Zi, f^ ti,tf '{zi) is essentially 
a continuous composition of the generator map v(z,t). 
Physically, it represents the effect of time evolution on 
the system, from the initial time ti to the final time 
tf, and is standard in classical field and semiclassical 
theories (T), HI, |45| . Note that the time-dependence of 
^(z, t) is analogous to classical dynamics under a time- 
dependent potential, thus requiring both superscript time 
parameters, ti and tf, rather than just the difference, 
At = (t f -t t ). 

If v(z,t) is presumed to be analytic, at least locally in 
the vicinity of some given trajectory, then f^ ti ' tf ^ l {zi) will 
also be locally analytic. Globally however, f^ ti,tf \zi) is 
usually not analytic even if v(z, t) itself is globally an- 
alytic, due to branch cuts that can arise in the func- 
tion (Zi). In most cases, e.g. when ^(x,t) has 
nodes, v(z,t) is not (globally) analytic, but "meromor- 
phic," meaning that it is globally analytic apart from 
a discrete set of simple poles [3g|. In this case also, 
f tyti ' tf \zi) is usually not meromorphic itself (though it 
must have simple poles), but also exhibits branch cuts. 
The branch cuts delineate regions of the Zi domain that 
get mapped discontinuously under f( tutf \zi) to different 
regions of the range. In the trajectory interpretation, two 
initially nearby trajectories, starting differentially close 
to each other at ti but on either side of a branch cut, 
wind up far apart at tf. The branch cuts themselves 
must therefore correspond to "separatrix"-type trajec- 
tory orbits (Fig.©. 

On the other hand, for trajectories sufficiently far from 
a separatrix or simple pole, f tyt " ts \z i ) may be regarded 
as locally analytic. Consequently, points z = Zi + dzi 
lying within a differentially small circle centered at Zi get 
mapped under (z) to another differentially small 

disk of points f^ ti,ts \zi + dzi) = Zf + dzf, as in Fig. [T] 
Due to local analyticity [36|, the ratio of differentials, i.e. 
dzf/dzi, is independent of the magnitude or direction of 
dzi. 

^ = /<*"*>'(*) (32) 
dzi 

Thus, the spatial derivative function, /(**•*/ ) (z^) specifies 
how the differential volume element, dz, is transformed 
via trajectory evolution over the specified time interval. 
For a given initial point, ti), we can, by varying tf = 



o 




position, *=Re[z] 



FIG. 1: Two closed complex trajectory orbits for a nonde- 
generate stationary state over half a period, At = T/2. All 
units are atomic units. Solid curve denotes main trajectory, 
starting at (z% = Xi,ti) and ending at [zf = Xf,tf). Dashed 
curve denotes neighboring trajectory, starting at (xi + dxi). 
The function f^'^^z) maps Zi to Zf, and the small differen- 
tial circular disk Zi + dzi on the right to the small differential 
circular disk Zf + dzf on the left. Right and left arrows de- 
note (directional) dxi and dx f intervals, respectively; note the 
change in sign. 

t, obtain a transformed dz for every point along the z(t) 
trajectory. In particular, 

^- = f^'(z i )=r^[z(t)}, (33) 

where it is understood that the trajectory, z(t), passes 
through (zi,ti). Note the new quantity, r tyZi ' ti ^[z{t)}, 
which — though equivalent to /(**'*) (z,j) — is introduced so 
as to be regarded as a function of the final trajectory 
point z(t) (at arbitrary time t), rather than the initial 
point Zi. Note also that A Zuti \z,j) = 1. 

For differentially small time increments, At = dt, it 
can easily be shown using standard statistical or hydro- 
dynamical arguments that dzf = [1 + v'(zi, ti)dt] dzi [Iol |. 
implying that 

^^M = ,'[z(t),<]r(-^)[z(t)] (34) 

The solution for A Zi ' ti%l [z{t)} is therefore found to be 

r (*.*0[ z ( t )] = ex p y v'[z{t'),t']dt' , (35) 

where the integration is along the trajectory. 

The general global solution to Eq. ([3"5]) above is not 
straightforward; in particular, it is not apparent how to 
obtain r^**'**' values for (z,t) points off of the z(t) tra- 
jectory. However, this is easily achieved for the specific 
case of nondegenerate stationary states, on which we will 
focus for the remainder of this paper. In this context, 
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v(z,t) = v(z) is independent of (final) time, and simi- 
larly, r^ 4 '*" = r Zi is independent of initial time, as ap- 
propriate for the resultant "time-independent potential" 
dynamics. A simple change of variables in Eq. (|35[) above 
then leads to 



r 2 *(z) 



v(zi) 



(36) 



Since v(z) = vim(z) is pure IMAG, for trajectories initi- 
ating on the real axis, Zj = Xi, the function r Xi (z) is pure 
REAL. Thus, whenever the trajectory z(t) recrosses the 
real axis, at (zf = xj,tf), the initial real- valued volume 
element dzi = dxi is transformed to a final volume ele- 
ment dzf = dxf that is also real- valued (Fig.[T]) — a non- 
trivial but essential prerequisite for any Eq. (|30p -tvpe re- 
lation, since p(xi, tj) and p[xf, tf) are also be real- valued. 
In particular, we obtain 



dxj 



dx 



f 



VlM^Xi) Vl M (Xf) 



(37) 



However, it is not necessary to restrict oneself to Xf val- 
ues that are connected to Xi via trajectories; in fact, the 
above relation is true for completely general values of Xi 
and x f . 

Substituting Eq. (J2(| into Eq. ([37]). and exploiting 
Eq. P5|) . we derive the following relation: 



dxj = 



P( x f) 



dx 



(38) 



Equation (|38p above is similar to, but clearly distinct 
from, the desired probability conservation condition, 
Eq. (|30|) . The latter is thus satisfied if and only if 
|p'(xj)| = |p'(a;/)|, which — though not true in general — 
is true for the special case of symmetric analytic velocity 
fields [v{—z) = —v(z)\ with a single stationary point at 
the origin. 

Note that for f At (zt) now depends only the time in- 
terval, A t . It is of interest to relate v(z) to f At (zi) via 
r Xi (zf). With z — Zi and / = f At (zi) — Zf, one obtains 



f =r^z f ) = V M 
dz v{z) 

w(f) — w{z) + const, 



(39) 
(40) 



where w(y) is the integral J dy/v(y). Clearly, w(y), and 
thus f At (z), will not in general be globally analytic, even 
if v(z) itself is. 

The role of stationary points, i.e. zq such that v(zq) — 
0, is important in trajectory dynamics. These occur 
where ip'(z) = 0, i.e. the non-node zeroes of p'(z) = 0. 
When ip(z) has nodes, it is easy to demonstrate that v(z) 
has a simple pole at each node, and is therefore mero- 
morphic at best. Conversely, node- free analytic wave- 
functions ip(z) can often lead to globally analytic velocity 
fields. In the vicinity of nodes, v(z) directs initially neigh- 
boring trajectories to very different final destinations, 



thus leading to separatrix-like trajectories and branch 
cuts in f At (z) (Fig. [2]). However, unlike separatrix orbits 
in classical phase space, v(z) does not approach zero as 
the node is approached — because the poles are not them- 
selves stationary points. Note that separatrix-like trajec- 
tories, and associated f At (z) branch cuts, can arise even 
when v(z) is globally analytic (Sec. IVI[) . 

In the neighborhood of stationary points, one can ap- 
ply the standard velocity field linearization method [45[ 
to determine the behavior of neighboring trajectories. In 
particular, v{z) w A\(z — zq), where A\ is a complex 
constant . If A\ is pure imaginary, then the neighbor- 
ing trajectories are circular and closed, with frequency 
u> = | Ai| and period T = 2ir/\Ai\. This is always the 
case for stationary points on the real axis, i.e. z$ = xq, 
which we now consider in greater detail. Note that 
due to reflection symmetry about the real axis, the first 
crossing of the real axis occurs at time T/2, at location 
Xf = xo — (xi~xo), for differentially small (xi—xo). Thus, 
in some small neighborhood of xq, the map f T ^ 2 (x) is 
real- valued, which in turn implies that f T ^ 2 (z) is REAL 
over the domain around xq for which f T ^ 2 (z) is analytic. 
This in turn implies that all trajectories within the do- 
main of analyticity have period T, even well away from 
the neighborhood where the velocity field is linear and 
the trajectory orbits circular. 

For stationary points off of the real axis, it may still 
be true that A\ is pure imaginary, leading to closed 
neighboring trajectories, and essentially "conservative" 
dynamics. However, A\ may also be complex-valued, re- 
sulting in aperiodic trajectories that spiral in or out, with 
the stationary point respectively an attractor or repeller 
(with the former akin to dissipative dynamics) . Such sta- 
tionary points zq come in complex conjugate pairs, with 
the Ai of one the complex conjugate of the other, thus 
implying that one point of the pair is an attractor and 
the other a repeller. It is possible for trajectories to flow 
directly from an attractor to its conjugate repeller, cross- 
ing the real axis in a region without closed orbits. 



VI. SPECIFIC EXAMPLES 

To demonstrate the range of dynamical behaviors 
available, both for analytic and non-analytic v(z), we 
consider a representative sampling of d = 1 systems, 
some of which have been considered in the previous liter- 
ature [HI, [I?], [H, US ■ In each case, the Hamiltonian 
is of the form 



h 2 d 2 



(41) 



with m — h — 1. All units may therefore be taken to be 
atomic units. Complex quantum trajectories for systems 
II- V are presented in Fig. [2] In every case, Eq. ([38]) has 
been confirmed numerically, for those trajectories that 
recross the real axis. The velocity field v(z) is obtained 
using Eq. (P0"]) ; stationary points are obtained by solving 
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(b) 



g 



■i 1 






1 / 


J-* / 
' i _!^s_ 





-1 1 

position. j:=Rs[;] 




-5 

position, jifRe[;] 



■S -2-10 1 
position, j-ReE^J 



FIG. 2: Complex quantum trajectories for a representative 
sampling of one-dimensional nondegenerate stationary states 
for Hamiltonians of the Eq. (|41[) form. All units are atomic 
units. Dashed lines indicated separatrix trajectories, leading 
to branch cuts in the map f T ^ 2 {z). See text for additional 
discussion. Specific states are as follows: (a) System II, Har- 
monic Oscillator First Excited State; (b) System III, Sym- 
metric Double-Peaked Ground State; (c) System IV, Morse 
Oscillator Ground State; (d) System V, Asymmetric Single- 
Peaked Ground State. 



System II: Harmonic Oscillator First Excited 
State 



potential: 
wavefunction: 
velocity: 
node: 

stationary point: 
period: 



V(z)=z 2 /2 
ip(z) = zexp(-z 2 /2) 
v(z) = i(z 2 - l)/z 
z = 
zo = ±l 

T = TT 



The first excited state of the harmonic oscillator has a 
node at the origin. As per the discussion at the start of 
Sec. El this implies that the velocity field v(z) is mero- 
morphic, with a simple pole at the z = node. Through 
this node, separatrices partition two sets of closed trajec- 
tories (period ir) around the two stationary points. Note 
that since these stationary points lie on the real axis, the 
nearby trajectories must be concentric circles (A\ is pure 
imaginary). There are also larger trajectory orbits, also 
closed (period 2ir), that surround both stationary points. 
Thus, trajectories on either side of the separatrices get 
mapped to very different locations, so for this example, 
f T / 2 (z) is not globally analytic, or even meromorphic. 



System III: Symmetric Double-Peaked Ground 
State 



v(zq) — 0; corresponding trajectory periods are obtained 
as per the end of Sec. [V] 



A. System I: Harmonic Oscillator Ground State 



potential: 
wavefunction: 
velocity: 
stationary point: 
period: 



V(z) =z 2 /2 

ip(z) = exp(-z 2 /2) 

v(z) = iz 

z = 

T = 2tt 



wavefunction: 
velocity: 
stationary point: 
period: 



ip(z) = cxp(-z 4 /2 + z 2 ) 
v(z) = 2iz(z 2 - 1) 
z = 0,±1 

T = 7T,7T/2 



There are no wavefunction nodes, and therefore no sim- 
ple poles in the velocity field v(z), nor in the trajectory 
map f T / 2 {z). Although v(z) is globally analytic, f T ^ 2 (z) 
is not. All trajectories except for the separatrices are 
closed, each surrounding exactly one of the three sta- 
tionary points. This is the first case considered for which 
there are separatrix trajectories when v(z) itself is glob- 
ally analytic. 



This system conforms to the special symmetric case, dis- 
cussed after Eq. [[58"]). for which |p'(xi)| = \p'(xf)\. Due 
to symmetry, Xf = —Xi corresponds to the first real- 
axis crossing at time T/2. Consequently, dxf = —dxi, 
p{xi) = p(xf), and weak probability conservation is sat- 
isfied. This is the only example considered for which 
Eq. ([30]) is correct. Also, rather uncharacteristically, 
f At {z) — e lAt z is globally analytic. The complex tra- 
jectory orbits are concentric circles centered at the origin 
that move counterclockwise. 



D. System IV: Morse Oscillator Ground State 



potential: V[z) = exp(-2\/2z/3) 

-2exp(-V2z/3) 
wavefunction: ip(z) = 6 18 1 / 4 

x exp^e-^ 2 / 3 - 5V2z/6) 
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velocity: v(z) = i(5 - 6e"^ /3 )/3V2 

stationary point: zq = 3 log(6/5)/\/2 + ik'iV2ir 

with k an integer 
period: T = 18tt/5 

The velocity field v(z) is globally analytic, but / T / 2 (z) 
is not. Also, v{z) is periodic in the imaginary direction, 
over a distance ?>^/2tt. Thus, the stationary point along 
the real axis, and surrounding "librational" closed tra- 
jectories, are duplicated at regular intervals away from 
the real axis. In addition, there is a family of open, "hin- 
dered rotational" trajectories, on the right side of the 
figure, that traverse all unit cells. 

E. System V: Asymmetric Single-Peaked Ground 
State 



wavefunction: "0( z ) — exp(— z 4 — z 2 /2 — z 3 ) 

velocity: v{z) = i(z + 3z 2 + 4z 3 ) 

stationary point: zq — 0, — 3/8 ± i\fl j% 

period: T = 2ir 

The velocity field v(z) is globally analytic, but / T / 2 (z) 
is not. The stationary point at z = is surrounded by 
closed trajectories with period 2tt, on the right side of the 
figure. In addition, there is a pair of complex-conjugate 
stationary points off of the real axis, such that the one 
in the upper half plane is a repeller, and the other an 
attractor. On the left side of the figure is a family of 
aperiodic, doubly spiraling trajectories that connect the 
repeller to the attractor. These trajectories cross the real 
axis only once. 

VII. CONCLUSIONS 

We conclude with a brief summary of what has been 
achieved here. First, it seems evident that any anal- 
ysis of quantum probability flux on complexified space 
requires generalized complex conjugation of the form of 
Eq. (|15[) . This results in analyticity of the requisite quan- 
tities such as probability density — but much more impor- 
tantly, leads to complexified flux relations [e.g., Eq. (f2"T)) ] 
that are physically relevant because they do not depend 
explicitly on the potential energy. On the other hand, 



the most straightforward complex generalization of the 
flux continuity relation to emerge from this work, i.e. 
Eq. (f2"T|) , corresponds to the v^E(z,t) rather than the 
v{z, t) velocity field. Even relaxing to the weaker con- 
dition of Eq. ([30]) is insufficient to achieve probability 
conservation for v(z,t). Yet even if this condition were 
somehow satisfied, it might not be very useful for most 
multidimensional applications, because few if any trajec- 
tories recross all of the real coordinate axes simultane- 
ously (even for separable systems, if the frequencies are 
incommensurate) . 

The above suggests that a different choice of com- 
plex velocity field might be more advantageous, though 
it is not clear at present how to construct such a field. 
One possible approach, inspired by coherent state initial 
value representations [l|, might be to regard z and z* as 
completely independent quantities, e.g. in formulating 
partial derivative expressions for the flux and its diver- 
gence. This might in principle require a pair of distinct 
complex trajectories, one on z-space and another on z*- 
space, which somehow get combined together to form the 
time-evolving probability density, p[z(t), z*(t),i\. Per- 
haps even two time coordinates would be involved. 

On the other hand, there is still some hope for the 
v{z, t) approach, as well. In particular, as per the dis- 
cussion following Eq. ([36]) . the fact that real dxi implies 
real dxf under v(z,t) dynamics might turn out to be 
quite useful. Also, in certain special cases, it may turn 
out that Eq. (|30p is approximately satisfied — well enough 
to enable calculations for large systems. In any event, 
it is hoped that this initial foray may enable subsequent 
developments towards achieving probability conservation 
for complex quantum trajectories, or at least, provides a 
useful framework for analysis of these important issues. 
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